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We propose octahedral group O^ as the family symmetry of neutrino-lepton sector. We find that 
O^ contains subgroups TJ^ (^Vi and Zi for realizing the bimaximal (BM) mixings, 623 = 0\2 = 45° 
and 013 = 0° , where Ui^ ® Z| and Z| serve as the residual symmetries of neutrinos and charged 
leptons, respectively. We present geometric interpretations of BM mixing in the octahedron, and 
construct natural geometrical breaking of Z| , leading to nontrivial deviations from the BM mixings. 
Our theory makes truly simple predictions of a relatively large reactor angle, ^13 ~ 45° — 612 = 
7.5° — 13.7° (3(t), the nearly maximal atmospheric angle and the approximate maximal Dirac CP 
violation. These agree well with the current neutrino data, and will be further probed by the on- 
going and upcoming oscillation experiments. 
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1. Introduction 

The observed three families of neutrinos, leptons and 
quarks in nature exhibit very distinctive and puzzling 
masses and mixing structures, posing great challenges 
to particle physics today. Unraveling the underlying 
family symmetry holds the best promise to overcome 
such challenges. For more than a decade, various oscil- 
lation experiments have discovered massiveness of neu- 
trinos and large atmosplieric and solar neutrino mixing 
angles {Q^z-! ^12)7 leaving the reactor mixing angle 6*^3 
much smaller and least determined. Inspecting the cur- 
rent oscillation data [l| shovifs that the atmospheric angle 
is around maximal, ^23 = 0(45°) , and the solar an- 
gle 0^2 is large, but has significant deviation from its 
maximal value, 45° — By2^ — 0(10°) . Strikingly, the 
recent T2K data [2] point to a nonzero O^^^ at 2.5(7 
level with central value Q^z = '^■'^° (11.0°) for normal 
(inverted) mass ordering, which is on the same order 
of the deviation 45° - 0^^ =■ 0(10°) . Subsequently, 
Minos [31 and Double Chooz [J| found further indica- 
tions of a nonzero 6^1, ■ The Daya Bay collaboration [5| 
newly announced a 5.2cr discovery of the nonzero O^i, , 
sin^ 26li3 = 0.092 ± 0.016(stat) ± 0.005(syst) , with cen- 
tral value 6']^3 = 8.8° and the 5cr allowed range. 
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Shortly afterwards, RENO [5| found nonzero dy^, at 4.9(7 
level, sin^ 26*13 = 0.113 ± 0.013(stat) ± 0.019(syst) , with 
central value 6*1, = 9.8° . 
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Both data are consistent with the earlier results of 
T2K 3, Minos [3, and Double Choozjl last year. 

From the available oscillation data [l|-[5|: we observe a 
striking pattern of neutrino mixing angles, 



where the second relation can be reexpressed in radian, 
6*13 = 0{\ — 9^2) = 0(0.2). Hence, we are strongly 
motivated to take the bimaximal (BM) mixing [6|[7[ as 
our leading order (LO) structure, 
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and then define the small deviations, (5^ = -|- — 6*33 , 
Sg = ^ — 9-^^2 7 9'iid (5^ = 6'i3 — , as the next-to-leading 
order (NLO) perturbation. 
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As shown before, in the lepton mass-eigenbasis, the 
neutrino sector exhibits a generic LO symmetry Iji^'^ (g) 

13) = (f , 0) and Zi 



Z| , where Z^'^ dictates 



determines 612 by its group parameter k = tan 9 12 (with 
k being a general real parameter of the 3-dimensional 
representation of Zj ) [8|. We find that the LO mixing 
pattern (j4|) is a specific realization of Zj^ 
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group parameter k = tan 6*12 = 1 ■ We also note that the 
previously much studied finite groups [7|[9|[l0|, such as 
Ai [ll| and 5*4 [1^, etc, do not contain Zj^ (g) Z|(A: = 1) 
as the subgroup. In this work, we demonstrate that the 
octahedral group O^ can unify Zj^ (E) li2{k = 1) and 
thus predict the LO structure ^. We further construct 
a natural geometrical breaking of the residual symmetries 
of O^ , leading to a quantitative prediction of the NLO 
structure ([5|) , as well as the maximal Dirac CP violation. 

2. Geometrical Formulation of BimELximal Mixing 

Most previous studies focused on the finite groups A4 
[111 ] and 54 [I4I, in order to realize the tri-bimaximal 
(TBM) mixing ansatz [7| which is a special realization 
of the generic product group Z2^ (g) Z| by choosing 
the group parameter of Z| to be, k = tan 6*12 = -y^ 
(6*12 — 35.3°), as shown in [8|. Another simple realization 
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FIG. 1. Geometric interpretations of A4, 5*4 and O^, as well 
as their relations. A4 only includes rotations of regular tetra- 
hedron and S4 is the full symmetry group of regular tetra- 
hedron, thus A4 C 5*4 . Moreover, S4 is a subgroup of O^ 
because the tetrahedron can be contained in a cube, which is 
the dual form of octahedron. The duality means a cube can 
be embedded into a octahedron with its 8 vertices being the 
face-centers of the octahedron, and vice versa. 
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whose 012 agrees with the current data even better than 
the TBM value. But any such choice of /c S Z| is inde- 
pendent of 7^2'^ (which dictates 9^^ = and 6. 
and thus docs not give the prediction of 61^ 
from these studies, the present work will start from the 
BM mixing Q as our LO structure, which corresponds 
to Z2^ (g) Z|(fc = 1) . We find the octahedral group O^ 
to be a minimal group that contains Zj^ (g) Z|(fc — 1) as 
subgroups. 

Since A4 and S4 arc much used in the literature [7| , we 
first analyze the relationships among A4, S4 and O^ . 
The A4 is the rotational group of a regular tetrahedron, 
and has 12 elements in total, which includes four 120° 
rotations plus their inverses, and three 180° rotations. 
S4 is the full symmetry of a regular tetrahedron and has 
24 elements, including the same rotations as A4 and ad- 
ditional reflection. These are interpreted by the two up- 
per plots in Fig.[T] The difference between A4 and ^4 
is the additional reflection in ^4, which corresponds to 
Z2^. Fig. [T] provides an intuitive geometrical explanation 
on why Zi^"^ appears accidentally in the A4 models. 

A similar situation happens for the BM mixing where 
the Z2^ appears accidentally when the BM mixing is ob- 
tained from an S4 model [9| . For this we also have a 
geometric interpretation, as presented by the lower plots 
of Fig.[T] The reason that an 6*4 model may obtain BM 
mixing is not due to its tetrahedral property, but because 
it is isomorphic to chiral octahedral group O, a subgroup 
of octahedral group O^ . The O,^ is the full symmetry 
group of a regular octahedron and isomorphic to O Zj 




FIG. 2. Assignments for the left-handed neutrino mass- 
eigenstates (|«^i") , \i^T) 1 l^^s")) ^-nd the left-handed charged 
lepton mass-eigenstates (l^™) , |^™) , |^™)) in the regular oc- 
tahedron, where Q, M, H, and G are the midpoints of the 
corresponding edges. 

or S'4 (g) Z2 , where Zj is a discrete inversion symmetry. 
The group O^ is well-known in crystal physics. 

Then, we present a geometrical formulation of the 
BM mixing ^ in the octahedral structure. Regarding 
the space in Fig. [5] as flavor space, we assign the left- 
handed neutrino mass-eigenstates (|j^™) , |i^™) , It'™)) 
and the left-handed charged lepton mass-eigenstates 
(1^™) , \i™) , |£™)) to certain vectors in the regular oc- 
tahedron with proper orientations. Thus, the transi- 
tion amplitude of \i^T) being an |i?™) -flavor neutrino 
is {iY^\i^T) ■ We can compute this via inner products of 
the corresponding geometrical vectors. 
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which is just the leptonic mixing matrix in the charged 
current, also called Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix IJ]. We readily recognize that (jG)) just 
results in the LO structure of BM mixing (j4]) . 

It is interesting to note that tetrahedron, cube and 
octahedron discussed above all belong to Platonic solids 
which are related to the classical elements of fire, earth 
and air, respectively, according to the ancient Greek 
philosopher Plato. The aesthetic beauty and symmetry 
of these objects have been appreciated by human since 
thousands of years ago. 

3. Application of O^ to Neutrinos and Leptons 

The full symmetry group of regular octahedron, O^, 
can be generated by two of its elements s and t with 



the definition relations, 



s^ = t^ = {stf = [tshsf = 1 
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The geometrical meanings of s and t are 90° rotoreflec- 
tioTi about the axis OE and 180° rotation about the axis 
OM, respectively (Fig. [3]). Here, by 90° rotoreflection, 
it means a 90° rotation around the axis OE and then 
making a reflection about the plane ABCD. 

The Oh has 48 elements which can be divided into 10 
conjugacy classes. Elements in the same class have simi- 
lar geometrical meanings. It has 10 irreducible represen- 
tations (IRs) and the sum of squares of their dimensions 



Dj equals the order of O/^ , 2, ^^ ^ 48 , which has the 

solution, D = {1,1,1,1,2,2,3,3,3,3}. This shows that 
Of^ has four 1-dimension (Id) IRs, denoted by 1, 1', 1, 1 ; 
two 2d IRs, denoted by 2, 2; and four 3d IRs, denoted by 
3, 3', 3, 3 . The IR 3 is the fundamental representation, 
under which the elements s and t take the forms, 
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Here for notational distinction, we use uppercase (low- 
ercase) letters for the representations (elements). The 
3d representations ([5]) can be directly inferred according 
to the geometrical meanings of s and t . In Fig.[3J we 
define OE as cc-axis, OA as y-axis, and OB as z-axis. 

The Of^ flavor symmetry breaking is expected to occur 
at a higher energy scale and leave the residual symme- 
tries in the low energy mass-matrices of neutrinos and 
leptons. For this work, we study the low energy effective 
theory and its phenomenological predictions, where the 
detailed high-scale dynamics of O^ breaking is irrelevant. 
In our effective theory, the flavor symmetry breaking is 
explicit (which is similar to the well-established explicit 
breakings of discrete P and CP symmetries in the SM) . 
Such breaking may arise from the spontaneous symme- 
try breaking via proper flavoron fields at high scales [7[ 
or due to certain unknown new mechanism, but the asso- 
ciated heavy states are all integrated out from the effec- 
tive theory and thus irrelevant to our low energy study. 
Since spontaneous breaking of discrete flavor symmetry 
still suffers other problems and complications including 
the well-known domain wall problem, we prefer explicit 
breaking for our effective theory. In our construction, 
the residual symmetry for leptons is Z4 , and is gener- 
ated by the 90° rotation r^ around the axis OE in Fig.[3l 
For neutrinos, the residual symmetry is "Zi^^ ® Zt , and 
is generated by the two reflections. 
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FIG. 3. Neutrino symmetry Zj^ (g) Z2 and charged lepton 
symmetry Z4 under octahedral group O^. The Zj^ corre- 
sponds to reflections through the plane PNQM (marked by 
hachured lines), and the Zl generates reflections through 
the green plane PAQC. Reflections through the three marked 
planes (hachured, green, blue) generate all symmetries in the 
neutrino sector. The Z4 corresponds to rotations by 90° 's 
around the axis FE. 

picted in Fig. [3] We can express them as below, 
r, = t{st)\ 



Zi: 

Zf: 

Zi: 



: StS'^tst , 

sts^ts . 



(9) 



Then, we derive their 3-dimensional representations. 
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In our construction, we assign the three generations 
of leptons to the triplet representation 3 . Then, we 
can write down the effective Lagrangian for lepton and 
neutrino mass-terms at the leading order, 

£0 = Ml^l^l + M2(^24 + ^3^3) + M3(^2^3 - ^3^) (11) 

+"hi(^i^i+ ^3^3) + "^22^^2^2 + 2mi3i/iJ/3 + h.c, 

where {1/1,1/21 '^3)'^ = ^ ^^^d {ii, i^, ^3)^ = i are left- 
handed neutrinos and leptons, and {i\, i\, ^3) = i'^ de- 
note the right-handed anti-leptons. The Lagrangian (fTTj) 
is invariant under the above residual symmetries, i.e., the 
Z4 for leptons, 



i -> R^i, 



"R 



t ' 



fl2) 



and Tl^ < 



' Z2 for Majorana neutrinos, 



(13) 



The mass matrices of leptons and neutrinos from ()lip 
can be diagonalized [l5| by the mapping matrices t/^g 
and Uuo, respectively. 
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where [/q is a diagonal Majorana-phase matrix absorb- 
ing possible phases from the neutrino mass-eigenvalues. 
Thus, we can deduce the PMNS mixing matrix. 
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UjM, 
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Uou;,: 



(15) 



where Uq exactly agrees to the BM mixing © (as de- 
rived earlier from the geometrical construction in Fig.[2|), 
and Uq = diag(l, 1, i)C/o is the Majorana-phase matrix. 
We also note that in the lepton-mass-eigenbasis, 
the symmetry transformation matrices of (jlOp become, 
(i?^, R'^^, R'^) == U^q{R^, R^^, RsW^q , with 

R'f, = diag(l, i, -i) , 
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Here, the transformation R' is just the conventional 
form of 2,2^ in the lepton-mass-eigenbasis and corre- 
sponds to the symmetry under the v and v^ exchange. 
As a final note of illustration, we give a type-I seesaw 
realization of the neutrino masses in (|lip . where we as- 
sign the Higgs doublet $ of the standard model (SM) as 
singlet of the Oj^ group, and the three heavy right-handed 
neutrinos [vi, i^f, v^) = v'^ as triplet, v'^ — > v'^Rj^^^ ^ under 
Zg^ (g) Z| transformation. Thus, we can write down the 
seesaw Lagrangian to replace the second row of (jlip , 

4s="hl('^li'l+'^3^3) + "^22^2^2 + "^13('^l'^3 + ^3^l) 

-A-h^ivlvl + vlvl)~ M22vlvl - '^M^zvlvl , (17) 

where the Dirac masses arise from the products of the 
corresponding Yukawa couplings and Higgs vacuum ex- 
pectation value ($) . Integrating out the heavy j/^ 's from 
(|17p . we derive the light neutrino seesaw- mass-terms in 
the effective Lagrangian (fTTj) with the coefficients, 
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At low energies our effective theory has the same particle 
content as the SM with light massive neutrinos. The 
heavy states associated with possible flavor breakings are 
integrated out around the seesaw scale, and thus cause 
no visible flavor changing effect at low energies. 

4. Geometrical Breaking and Physical Predictions 

Sec. 3 demonstrated how the LO BM mixing (|4]) is pre- 
dicted by the residual symmetries ll^ ® Z| and Z4 of 
the octahedral group O^ . In this section, we will con- 
struct a natural geometrical breaking of Z4 in the lepton 
sector, and derive physical predictions. 

Inspecting the residual symmetries Tl^ ® U^ ^'^d ^4 j 
we see that there are two fundamental ways to introduce 
the breaking: one is to break the Z4 in the lepton sector, 



and another is to break TL^ 



) Z2 in the neutrino sector 



Mfi-Mf3 



(as in [8|[13|). In either way of the breaking, it is gen- 
erally recognized that charged leptons always have the 
maximal residual symmetries U [\) ®U (\) ®IJ iX) 3 Z4, 
and neutrinos always have the maximal residual sym- 
metries Z2 (81 Z2 (8) Z2 D Z2 (8) Z2 , even though after 
any possible breaking the Z4 will differ from the origi- 
nal Z^ and Zg ® Z^ is no longer the Z^^ ® Z| . The 
breaking of the original Z| or TJi^ ® Zj always accom- 
panies the emergence of a new Z4 or Z2 (8 Z2 . Hence, 
adding a residual symmetry breaking term (5£ to the 
LO Lagrangian £g just means changing the correspond- 
ing residual symmetry from Z4 to another Z4 , or from 
Z2^ 8) Z2 to another Zj 81 Zj . It usually turns out that 
the emergent new residual symmetries in the new La- 
grangian £g -|- i5£ , such as Z4 or Z2 81 Z2 , may be no 
longer physically useful, but the modiflcation of the orig- 
inal symmetry predictions based on the LO Lagrangian 
£0 due to the residual symmetry breaking term bH, docs 
lead to physically testable corrections. 

For our study, we find that the simplest way comes 
from a natural geometrical breaking of Z4 in lepton sec- 
tor. This geometrical breaking is constructed by simply 
rotating the axis OE of Z4 to another axis OE' asso- 
ciated with a new Z4 , where we define the small rota- 
tion angle ZEOE' = -\/2e with a characteristic breaking 
parameter e . We illustrate this geometrical breaking in 
Fig.m where the axis OE rotates to OE' within the plane 
GIJH and G (H) is the middle point of the line AD (BC). 
Thus we can derive the Z4 breaking Lagrangian, 

bL = -e[(/ii-A^2-/i3)^i^2 + (/i2-Ml-M3)^l^3 (19) 

+ (Mi-M2 + Ai3)^2^i + (Ai2-Aii + ^^3)^3^1] + h-c-, 

where /^i 2 3 have appeared in the LO Lagrangian (jlip . 
We find that the small breaking parameter e indeed gen- 
erates the NLO deviation via e ~ 45° - Q^^ =■ 0(0.2) . 
We have verified that Lq -\- 5C is invariant under trans- 
formations of i?« € Z4 . Geometrically, i?„ is a rotation 
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FIG. 4. Geometric breaking of Z4 is realized by pulling the 
Z| axis OE to OE', where the point E' is in the line IJ. 

around the axis OE' by 90°, and we deduce up to 0(e), 




FIG. 5. Prediction of 613 as a function of 45° — ^12 , with 
300 samples (marked by blue +). The yellow (pink) region 
corresponds to the 3a (5cr) new limit of Daya Bay [3. The 3a 
(5a") new limit of RENO 0] is given by the blue (red) dashed 
lines. The black dashed box depicts the allowed 3cr ranges by 
the global fit Q. 
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range of solar angle, 31.3° < 9i2 < 37.5° , with a central 



value 9 



12 



34.5° . With this, our first relation of 



predicts the 3a range of the reactor angle. 



Combining the breaking term ([T9|) with the LO La- 
grangian (jlip . we derive the full PMNS mixing matrix 
up to the NLO, 
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where C/' = t/g ^^ given below (|15|) since the symmetry 
breaking Lagrangian (|19p does not affect the neutrino 
mass matrix. It is clear that V differs from the LO result 
([T5|) by 0(e) corrections. We can also break the Z| by 
rotating OE in the direction perpendicular to the plane 
GIJH, but the difference can be absorbed by flipping the 
sign of Dirac CP phase. Finally, from (pTj) . we derive the 
following solution, up to 0(e) corrections. 
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where as we defined before, the three NLO parameters 
(<5„<5„JJ ^ (f "^23, f-^12,^13), and S^ is the 
Dirac CP phase angle. From (|22|) . we infer the truly 
simple and beautiful predictions below. 
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The first relation above strikingly predicts the crucial 
reactor angle 9^^ in terms of the precisely measured solar 
angle 0^2 ^ without explicit dependence on the breaking 
parameter e itself. The current global fit [l| gives the 3a 
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(24) 



This is rather distinct and agrees well with the new data 
from Daya Bay and RENO reactor experiments [5|. 

We further plot the predicted 9^^ as a function of 45°— 
0]^2_iii Fig-El where we have input 9i2 from the global 
fit [l| by generating 300 random samples with Gaussian 
distribution. Fig.[S] also shows the 3a constraint by the 
global fit [l| on the allowed parameter space, as marked 
by black dashed box. We further present the 3a and 5a 
limits of 9^^ from the new Daya Bay data[5[, as given 
by the horizontal yellow and pink bands, respectively. 
The 3a {5a) limits of RENO data 3 are shown by the 
horizontal blue (red) dashed lines. Our prediction on 9i^ 
(marked by blue +) lies in the upper half of the allowed 
regions by Daya Bay and RENO data, and will be further 
tested by the upcoming data of Daya Bay [5| , T2K [2] , 
Double Chooz [i| and RENO [l3] experiments. 

We also note that under the seesaw generation of light 
neutrino masses via (|17p - (|18p . it is natural to have the Oh 
breaking and Z4 breaking realized at or above the seesaw 
scale. Then, we can further include the small effects of 
renormalization group (RG) running for mixing angles 
([251). We find that the RG running effects are possible 
to reduce 6*33 and 9i2 + ^13 by about 2° at low energy. 
Thus, we have 9^3 ~ 43° and 9^^ ~ 43° - 6*12 , which 
will agree better with the recent global fit [l| . 

Finally, we further compute the Jarlskog invariant J 
[iq as a measure of leptonic CP violations. With 
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FIG. 6. Prediction of Jarlskog invariant J as a function of 
45°— 6i2 , with 300 samples (marked by red +). The yeUow 
band represents the 3a allowed ranges of 45° — 612 from the 
global fit [I[. 



we derive J in terms of the solar mixing angle 9i2 alone, 
J = -g- sin 2^23 sin 20]^2 sin 20]^3 cos ^]^3 sin (5j3 

~±j^ sin 4012- (25) 

From the global fit [1[, we derive the 3a range of the 
predicted Jarlskog invariant, 

\J\ = 0.031-0.050. (26) 

With Eq. ([25]) and inputting the global fit[l[ of 0^2 j we 
plot om- prediction of J as a function of 45° — 612 in 
Fig. [51 where the yellow region shows the 3cr limits on 
45°— 6*5^2 from the global fit. Our distinctive prediction 
of the nearly maximal Dirac CP violation will be further 
probed by the upcoming long baseline oscillation exper- 
iments. 



5. Conclusions 

Understanding flavor symmetries behind three fami- 
lies of the standard model fermions in nature holds the 
best promise to unravel their mysterious masses and mix- 
ing structures. In this work, we proposed octahedral 
group O^ as the flavor symmetry of neutrino-lepton sec- 
tor, whose residual symmetries Zj^ (g) Z^ and Z4 pre- 
scribe the neutrinos and charged leptons, respectively. 

The current neutrino data strongly motivate us to con- 
sider the bimaximal mixing (BM) of (JH) as a good LO 
structure and the deviations from the BM mixing in 
Eq. ([5]) can nicely serve as the NLO perturbation. In 
Sec. 2-3, we presented the geometric interpretations of 
Of^ and its residual symmetries (Fig. [1] and Fig. [3]). We 
found that the unique octahedral group O^ naturally in- 
cludes the subgroups Zj^ (g) Z| and Z4 for generating 



breaking of Z4 as in Fig. |4] and Eq. (fT9|) . The resultant 
PMNS mixing matrix V was derived in Eq. pT|) . Our 
theory makes truly simple new predictions of a relatively 
large reactor angle, 6*13 ~ 45° - 6*12 = 7.5° - 13.7° (3cr), 
the nearly maximal atmospheric angle, and the approx- 
imate maximal Dirac CP violation. These are shown in 
Eq. (|23p and presented in Fig.[5l Especially, our new pre- 
diction of the reactor angle 61^ agrees well with the new 
Daya Bay discovery dJ) and the RENO limit ^. Finally, 
we predicted Jarlskog invariant J in terms of a single 
mixing angle 612 via Eq. (|25|) . With this we derived its 
3cr range, \J\ ^ 0.031 - 0.050, as in Eq. ^. The cor- 
relation between Jarlskog invariant J and the deviation 
45°— 6*23 is depicted in Fig.[Bl All our new predictions 
agree well with the current neutrino data, and will be 
further tested by the on-going and upcoming oscillation 
experiments. 
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